Interaction of Ising-Bloch fronts with Dirichlet boundaries 
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We study the Ising-Bloch bifurcation in two systems, the Complex Ginzburg Landau equation 
(CGLE) and a FitzHugh Nagumo (FN) model in the presence of spatial inhomogeneity introduced 
by Dirichlet boundary conditions. It is seen that the interaction of fronts with boundaries is similar 
in both systems, establishing the generality of the Ising-Bloch bifurcation. We derive reduced 
dynamical equations for the FN model that explain front dynamics close to the boundary. We find 
that front dynamics in a highly non-adiabatic (slow front) limit is controlled by fixed points of the 
reduced dynamical equations, that occur close to the boundary. 

PACS numbers: 82.40.-s,05.70.Ln 



I. INTRODUCTION 



In spatially extended reaction-diffusion systems far 
from equilibrium, the interplay of the diffusion and re- 
action processes is frequently associated with the forma- 
tion of spatial or temporal patterns in the concentration 
fields pi 0, IE 0> IE El • O ne such example is a front-like 
structure connecting two different homogeneous steady 
states. In a bistable system, where both steady states are 
stable against small perturbations, these fronts can un- 
dergo bifurcations, known as nonequilibrium Ising-Bloch 
bifurcations, where a stationary Ising front exchanges 
stability with a pair of counterpropagating Bloch fronts. 
This bifurcation has been observed in several chemical 
reactions 0, |E 0] and also in liquid crystals 0,11,13 sub- 
ject to an external time-dependent perturbation. 

Two models of this Ising-Bloch bifurcation have been 
extensively studied in this context. One is the para- 
metrically forced Complex Ginzburg Landau equation 
(CGLE) 0. This system describes nematic liquid crys- 
tals subjected to a rotating magnetic field and a high 
frequency electrical field Q. The CGLE is often used 
in spatially extended systems to describe the dynam- 
ics close to an oscillatory instability (Hopf bifurca- 
tion]. The other is a FitzHugh Nagumo (FN) model 
H Hi [H [H EU^which qualitatively models various 
chemical reactions E3- E3- Front solutions in this 
model have been extensively investigated, specially when 
translational invariance is broken by the presence of spa- 
tial inhomogeneities, which is often the case in realistic 
experimental situations. This includes the ways in which 
one Bloch front can be perturbatively chan ged to the 
other (leading to front reversal) [ToL ITTl IL% Il3j . Par- 
ticularly, one such scenario for front reversals and other 
exotic nonuniform (variable velocity) front motion like 
breathing, involves the breaking of translational invari- 
ance by zero flux boundary conditions £|. This nonuni- 
form motion of fronts is explained by the presence of uni- 
form velocity front solutions (nullclines of the FN model 



partial differential equations) to which faster or slower 
moving fronts relax adiabatically. 

In this paper we examine nonuniform front motion in 
the case when translational invariance is broken by im- 
posing fixed chemical concentrations at the boundary of 
a reactor (Dirichlet boundary conditions). An adiabatic 
description of nonuniform front motion along the lines 
of |J] is inadequate. It relies on nullclines and fails to 
explain the jumps from one nullcline to another or the 
influence of fixed points. Therefore, we employ a dynam- 
ical approach which satisfactorily accounts for jumps and 
fixed point influence. We also establish the generic na- 
ture of nonuniform front motion in the presence of spatial 
inhomogeneities by studying it in two distinct systems, 
the FN and CGL models. 

The interaction of traveling fronts with boundaries 
for Dirichlet boundary conditions shows several new fea- 
tures. We see a transition from front reversal to trapping 
of an incoming Bloch front on its approach to a bound- 
ary as a function of boundary values. We also find that 
trapped fronts and reversed fronts can coexist for certain 
boundary conditions. Finally, we derive reduced dynam- 
ical equations that explain the features mentioned above. 

In section II of this paper we review the relevant details 
of the CGL and FN models. Section III presents our 
numerical study of the two models. Section IV contains 
an analytic study of front interactions with boundaries 
for the FN model. Our conclusions are listed in Section 
V. 



II. THE MODELS 

The parametrically forced CGLE is the generic model 
describing the slow phase and amplitude modulations 
of a spatially distributed assembly of coupled oscillators 
near its Hopf bifurcation 0]. This assembly of auto- 
oscillators is parametrically forced at twice its natural 
frequency and can be written as, 
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The complex field A contains the amplitude and phase 
of the coupled oscillators, fi measures the distance from 
the oscillatory instability threshold, v is the detuning of 
the forcing term frequency from exactly twice the Hopf 
frequency, a and (3 are real control parameters, and 7 > 
is the forcing amplitude at twice the natural frequency. 
The right hand side of Eq. cannot be written in a 
variational form if v, a, and (3 are nonzero. 

The parameter k represents parametric forcing at the 
natural frequency of the system. If k = 0, Eq. (|TJ has a 
the parity symmetry A — ► —A, and the noncquilibrium 
Ising-Bloch bifurcation is a symmetric pitchfork 0, 0] 
with the front velocity as the order parameter. Zero 
velocity Ising walls lose stability to counterpropagating 
Bloch walls, as the bifurcation parameter 7 crosses its 
critical value. The pitchfork unfolds into a saddle node 
bifurcation for a nonzero k where along with the stable 
Ising wall, a stable and unstable pair of Bloch walls ap- 
pear at the bifurcation. 

The CGLE Eq. has trivial homogeneous solutions 
(A(x) = iy/n + 7 + Oi) in the variational case (a — (3 — 
v = 0) making it bistable. Solutions connecting these 
trivial states are, 



A(x) = sfji + 7 tanh 



'M + 7 



•Oi, 



which is the Ising wall solution characterized by a zero 
imaginary part, and 



A(x) = + 7 tanh (y/2jx) + io\J [i — 37sech(y / 27x), 

the two Bloch wall solutions distinguished by their re- 
spective chirality a = ±1. In the case when nonva- 
riational parameters (a, (3, v) are non-zero, and 7 2 > 

solutions, 
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-(3 2 ), Eq. (JTJ is still bistable and has trivial 
i?exp(i0) 

M + /3j/+ [ 7 2 (1+ ^ ) _ (ly _ /3M) 2 ] l/2 



(1 + /3 2 ) 



cos(20) 
sin(20) 



(- M + £ 2 )/ 7 , 
{v-(3R 2 )h. 



(2) 



Nonvariational Ising and Bloch wall solutions are qualita- 
tively different from their variational counterparts. Non- 
variational Bloch walls move as a result of the breaking of 
chiral symmetry Q, whereas variational Bloch walls are 
stationary. Chirality breaking is unlike other mechanisms 
of front motion, where a glo bally stable state invades an 
unstable or metastable state [l!i 13] • In the perturbative 
limit when (a, [3, v) are small, the velocity of Bloch fronts 
can be written as 
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2(3/i- 7 ) 



\-v + ^+( a - pyj\ , (3) 



where \ — ±\/3(7 c 7) • The Bloch wall velocity c 
is proportional to V7c ~ 7 as expected for a pitchfork 
bifurcation. 



Now we look at the FN model, which also shows 
a front bifurcation. This is a simple two component 
model and has been thoroughly analysed in Ref. 0, 
IH 13 13 EH in the context of this bifurcation. It 
has been widely used to model patterns in reactions 
like the Belousov-Zhabotinsky (BZ) reaction [3 13 13 
13 H3- Ferrocyanide-Iodate-Sulfite (FIS) reaction ~Q 
and Chlorite-Iodide-Malonic-Acid reaction(CIMA) [lj, 
llH HtH . The two component reaction-diffusion system, 
with v(x,t) impeding the production of u(x,t), is given 
by 



du 

~dt 

dv 

at 



e 1 (u — u 3 — v) + S 1 u x 
(u - aiv - a a ) + v xx . 



(4) 



The parameters e and 6 differentiate the time scales and 
space scales of the two fields respectively. The parame- 
ters a\ and a characterize renormalized local reaction 
parameters, possibly after an adiabatic elimination of 
faster reacting species. Equations 10} are in general non- 
variational except for certain specific parameter values. 
The parameter a is analogous to the parameter k in the 
CGL equation, and it controls the symmetry of the front 
bifurcation. If a — 0, the FN model undergoes a sym- 
metric front bifurcation represented by a pitchfork. For a 
nonzero a a the pitchfork unfolds into a saddle node as in 
the CGLE. A notable difference between the two systems 
is the presence of the parameter 8 in the FN model. This 
parameter affects the relative spatial extent of the fronts 
connecting the trivial homogeneous solutions of Eq. 
Thus, by choosing a suitable e/5 ratio, the connecting 
fronts of one of the fields can be made very sharp com- 
pared to the other. This is not possible in the CGLE, 
where fronts for both the real and imaginary parts have 
the same spatial extent. 

Trivial homogeneous solutions to the FN model are 
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CL\ — 1 
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u c 



(5) 



Similar to the CGLE, the FN model has Ising and Bloch 
fronts as its solutions, that bifurcate in a pitchfork. The 
bifurcation parameter in the FN model may be chosen 
to be 77 = y/eS. The critical value of this parameter 
is Vc = (3/VT2)2)(ai + l/2) 3/2 [3] . The Bloch wall 
velocities given by 



6V2 



^01 + 1/2) 



(Vc - v) 



(6) 



are proportional to \Jr\ c — i]> the deviation of the bifur- 
cation parameter from its critical value, as expected for 
a pitchfork bifurcation. 

As discussed above, both these models have common 
features associated with the front bifurcation. This forms 
the basis of their comparative study in the forthcoming 
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sections. Front dynamics in both these systems can be 
represented by the system of equations, 



x = c 

c = Oc - p)c - gc 3 



(7) 



These equations employ the pitchfork bifurcation nor- 
mal form with velocity as the order parameter, coupled 
with the trivial observation that the velocity is the rate 
of change of position. The front velocity c and position 
x, therefore constitute two degrees of freedom, that are 
sufficient to describe front dynamics close to the front 
bifurcation. The bifurcation parameter p is denoted by 7 
for the CGL equation and rj c for the FN model. The two 
independent variables c and x in Eq. J2J, which are obvi- 
ously uncoupled, represent dynamics where translational 
invariance is present, and the solutions are independent 
of a choice of spatial origin. Scenarios can be envisioned 
where this translational invariance is broken. Examples 
are the imposition of different boundary conditions Q 
or spatial inhomogeneities introduced externally via an 
advective field [T^]- To account for broken translational 
invariance in these scenarios, one has to modify Eq. JJJ. 
This would lead to the coupling of the two degrees of 
freedom x and c, the details of which would depend on 
the scenario considered. In this paper, we introduce a 
spatial inhomogeneity in the form of Dirichlet boundary 
conditions, examining its effect on incoming Bloch fronts 
in both models, and drawing parallels. The next section 
constitutes our numerical study. In section IV, we derive 
and analyse the way in which Dirichlet boundary condi- 
tions couple x and c for the FN model. Such an analysis 
is not possible for the CGL equation since it is impos- 
sible to represent the front position by a single point. 
Therefore, we rely exclusively on numerical simulations 
for that equation. 
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FIG. 1: (a) Typical traveling front for the CGLE (b) Traveling 
front in the FN model 



III. NUMERICAL RESULTS 

In this section we lay out the numerical details of the 
study of the CGL and FN models. The simulations in 
the two models are carried out in regimes where analyt- 
ical calculations performed by reducing front dynamics 
to a fewer degrees of freedom are not possible. This is 
the regime where Bloch fronts have high velocities (far 
beyond the front bifurcation threshold), and the fields 
forming the front core, have similar spatial scales. 

We solve both the CGL and FN system of reaction 
diffusion equations using an implicit Crank-Nicholson 
scheme. Dirichlet boundary conditions are used at both 
ends of the domain, which is typically composed of 400 
grid points with a time step size of 0.01. The boundary 
values at one end are fixed at one of the homogeneous 
solutions of the Eq. JIJ and |@J| . At the other end we are 
free to vary the boundary condition. In our numerical 
simulations, we keep the domain large compared to the 
characteristic spatial extent of the front, so that the in- 
fluence of the boundary is only felt when the front is close 



enough to it. We verified that the grid and time steps 
were small enough to ensure that the numerical solution 
converged. 

By a suitable choice of initial conditions a Bloch front 
or its counterpropagating partner can be generated. The 
symmetry of the bifurcation ensures that they have the 
same speed as long as they are not close to a boundary. A 
typical front for the paramctrically forced CGL equation 
and the FN system of equations is shown in Fig. ^ Bloch 
fronts for the CGL model show a characteristic chirality 
broken structure at their core which is essential for 
their propagation. Similar structure considerations apply 
to fronts in the FN model [2l| . 

In our simulations we focus on the interaction of incom- 
ing Bloch walls with boundaries, where Dirichlet bound- 
ary conditions are imposed on the two fields, He A, ImA 
in the CGL equation and u, v in the FN model. Par- 
ticularly we look at how the front cores are perturbed 
by the boundary for a whole range of boundary condi- 
tions. Fronts coming into the boundary from infinity in 
both models either rebound or get trapped depending on 
the Dirichlet boundary values. A front that traps loses 
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FIG. 2: (a) A typical Nontrivial stationary solution to which 
trapped fronts evolve into for the CGLE. The spatially ho- 
mogeneous solutions are connected to the Dirichlet boundary 
values, (b) Nontrivial stationary solution for the FN model. 



its core structure and evolves into the nearest available 
stable (attracting) configuration of the fields, which in 
this case is the nontrivial steady state solution of Eq. 
or Eq. (J2J for that particular boundary condition. Fig- 
ure [3 shows a typical nontrivial steady state solution for 
the CGLE and FN model. This solution is comprised of 
a spatially homogeneous part, given by Eq. for the 
CGLE or Eq. JSJ) for the FN model, and an inhomo- 
geneous part that connects the spatially homogeneous 
solutions to the Dirichlet boundary value. Rebounding 
phenomena close to the boundary is characterized by the 
core of an incoming front flipping into the core of its 
counterpropagating partner, resulting in the front mov- 
ing away. 

Our observations are plotted in the plane of boundary 
values, revealing a curve separating regions of bouncing 
and trapped fronts for both CGL and FN systems. Fig- 
ure |3 shows the transition curve for a = —0.1, 7 = 0.31, 
(3 = —0.15, [i = TO and v = 0.1 for the parametrically 
forced CGL equation. A similar transition curve for the 
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FIG. 3: Transition Curve from a region of trapped Bloch 
fronts to bouncing ones for the CGLE. ReA boundary value 
on the X axis, ImA boundary value on the Y axis 
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FIG. 4: Transition Curve from a region of trapped Bloch 
fronts to bouncing ones for the FN model. U boundary value 
on the X axis, V boundary value on the Y axis. 



FN model, with a\ = 2.0, 8 = 0.14, e = 0.05 is shown in 
Fig. As one closes in on the curve from the trapping 
(bouncing) region, the fronts take longer to get trapped 
(bounce). This slowing down in the dynamics close to 
the transition curve is indicative of critical behavior, usu- 
ally associated with eigenvalues of fixed points that ap- 
proach zero when a parameter (the Dirichlet boundary 
value here) is varied. An analytical explanation of the 
slowing down is given in the next section. 

Incoming Bloch fronts, in both the CGL and FN mod- 
els, evolve into nontrivial stationary solutions to Eq. l|T]l 
and Eq. (@J) respectively when trapped at the boundary. 
These nontrivial solutions are linearly stable by virtue of 
the Bloch fronts evolving into them. It remains to be 
ascertained whether these solutions remain stable when 
Dirichlet boundary values that lead to a bounce are im- 
posed at the boundary. Hypothetically, one could as- 
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sociate the loss of stability of these nontrivial solutions 
with the transition from trapping to bouncing of incom- 
ing fronts. We test this hypothesis by carrying out a 
numerical linear stability analysis of the nontrivial sta- 
tionary solutions ipo- Eq. is linearized to, 



^ = (1 + ia)V 2 Sip + [/i + %v- 2(1 + i/3)\ip \ 2 ]5il) 
at 



+h - (i + w 2 a ]sr 



(8) 



for the perturbation Si/) = at the boundaries. A similar 
linearization is done for the FN system. It is found that 
the eigenvalue spectrum for both systems has negative 
real parts in the trapping region, as expected. More- 
over, these real parts remain negative when we evaluate 
the stability of the nontrivial stationary solutions in the 
bouncing region. Even though nontrivial solutions are 
stable in the bouncing region, incoming fronts do not 
evolve into them. Consequently, the critical behavior is 
not governed by the loss of stability of these solutions. 

We also observe that close to the transition curve 
bouncing and trapped fronts can coexist. Instead of 
coming in from infinity and rebounding, a Bloch front 
created close to the boundary may get trapped even if 
Dirichlet boundary values that produce a bounce are im- 
posed. To demonstrate this, we choose a boundary value 
inside the trapping region, close to the transition curve. 
A Bloch wall is launched from infinity towards the bound- 
ary. As the wall approaches the boundary, the simulation 
is stopped. The field configuration is saved and used as 
the initial condition for the next simulation run. In this 
new run we make a small change in the Dirichlet bound- 
ary value and move across the transition curve into the 
bouncing region. If the front core in the saved configu- 
ration is close enough to the boundary, the front will get 
trapped, even if Dirichlet boundary values that produce 
a bounce are imposed. An analytical explanation of co- 
existence for the FN model is given in the next section, 
which gives insights into the coexistence behavior in both 
models. 

Summarizing, in both systems we have regions of 
trapped incoming Bloch fronts and bouncing Bloch fronts 
in the plane of boundary values. Critical slowing down 
of the front dynamics in proximity to a boundary is ob- 
served in both systems close to the transition curve. Also, 
the nontrivial stationary solutions in both systems re- 
main stable across the transition curve, implying that 
they are not responsible for the critical behavior we see. 
Both systems exhibit the coexistence of bouncing and 
trapped solutions. 

In the next section, we explain these numerical obser- 
vations by deriving the mechanism that shows how x and 
c are coupled for Dirichlet boundary conditions in the FN 
model. 



IV. ANALYSIS 

The goal of this section is to explain front interaction 
with boundaries, in terms of coupling of the evolution 
equations for the two degrees of freedom, front velocity c 
and the front position x in Eq. J7J. It is shown that the 
coupling is the result of the spatial inhomogeneity sensed 
by the front as it approaches the boundary. 

We restrict ourselves to the regime where e/8 is small 
and rj — ^/TS is finite. This restriction leads to a very 
sharp spatial variation of u(x, f) field at the core of the 
front. The slowly varying v(x,t) field can be considered 
constant in this core region. Hence the v(x,t) field is 
represented by a single value Vf, at the point where the 
u(x,i) field has zero value. We also restrict ourselves 
to small front velocities c, which can be done, by either 
making rj c — rj small or choosing a large a± . 

The restrictions described above lead to a drastic re- 
duction in the number of degrees of freedom used to de- 
scribe the front. As opposed to a front description based 
on the whole model Eq. (0J, a sharp front (fast spatial 
variation of u at the core) can be thought of as a point 
particle with a definite position and velocity. The slowly 
varying v variable can be thought of as a field associated 
with this particle that allows it to sense the boundary. A 
small velocity has a few simple implications. An addition 
of two kinds of perturbations, one that changes Vf (veloc- 
ity) slightly and the other that produces a local distortion 
in the v(x,t) field far from the front position, to a slow 
moving and uniformly translating front, is followed by the 
disappearance of the distortion and the relaxation of the 
front back to constant velocity. The time scale on which 
the distortion vanishes is much faster than the scale on 
which Vf relaxes back to its original value. In essence, 
we have two time scales, the slower one associated with 
nonsteady front motion. 

We now employ the restrictions mentioned above to 
obtain a reduced description of Eq. J2J . We solve Eq. Q 
with a Dirichlet boundary condition v = Vb at the left 
boundary. On the right boundary, which is at infinity, we 
have v = —q~ 2 , and q 2 = a\ + 1/2. Following Ref. [l 
we have the following system of equations, 



(9) 



and 



v t +q v-V r 



Vt + q V — v r 



r]V2 
3 



vV2 

v(—x,t) = Vb , v(oo,t) = —q~ 2 



v(0, t)v r + 1 r<0 
v(0,t)v r -l r>0 



(10) 



Eq. © implies that the velocity is proportional to the 
value of the v(x, t) field at the sharp interface formed by 
the u(x, t) field. Eq. (|10fl represents the evolution of the 
v(x, t) field in a frame of reference which moves with the 
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front. The variable r is the spatial coordinate from the 
front position and r — —x is the distance of the boundary 
to the left of the front. 

We solve Eq. ijlO|) perturbatively. The starting point 
of the perturbative expansion is to find a stationary Ising 
wall solution. Therefore, setting the time derivatives in 
Eq. IjlOl) to zero and looking for an Ising wall solution, 
we get, 



q 2 v 



1 = 



q 2 v - 1 = 



r < 
r > 



(11) 

with v(0+) = v(0-) = and v(oo) = —q~ 2 . This ensures 
that that the Ising wall has Vf = 0. The solution to 
Eq. dU is, 



,(°) - 



,(0) 



= -q-\e qr -l) r<0, 
1) r>0. 



-2( p -qr 



q~ 2 (e 



Hence for the Ising wall at r = —x, we have, v^°'(—x) = 



(1 



C )/Y 



We look for traveling Bloch wall solutions as a pertur- 
bation to this uniquely defined Ising wall. Since the Bloch 
walls have a Dirichlet boundary condition v(—x) = Vb, 
the perturbative correction to the Ising wall should have 
a boundary value v c = V), — (1 — e~ qx )/q 2 which changes 
as the front moves. 

Let v be the perturbation. Then, 



,(o) 



(12) 



v is expanded in powers of c, the small perturbation pa- 
rameter. Since the front bifurcation is a pitchfork, 7/ is 
expanded in powers of c 2 . 



oc 

v(r,t,T) = J2 c(n)v(n) (r,t,T) 

n=l 

T] = i] c (x) - c 2 rji(x) + c 4 rj 2 (x). 



(13) 



T = c 2 t is the slower time scale responsible for nonsteady 
front motion. The coefficients of powers of c in the ex- 
pansion of r\ are functions of x to incorporate the broken 
translational invariance. Using Eq. I|12(l and Eq. (f 131) in 
Eq. (|1C)|> . one obtains 



,(n) 



2 „(») 



q v 



1,2,3.., (14) 



with 



.,(2) 



,(3) 



V2v, 
3 

..(1) . 



.J 1 ) v (o) 



J 1 ) y(l) 

u \r=0 u r 



3i?i (i) 

— v. 



(2) 



V \r=0 V r 



.,(0) 



J 1 ) v m +w (2) „ 



(2) 4_ ,.(3) ,.(0) 



.(15) 



We use Green's functions to solve the system of equations 
above. The general solution, given that we have found 
an appropriate Green's function G(r, t\r', t') is, 



v (n) {r,t) = j I G(r,t\r',t')p n (r',t')dt'dr' 

rt 



U Js' 
I 



d'v^lr' t') 
G(r,t\r',t') { ' ' 



d'G(r,t\r',t') (n) 

dn [r,t) 



dn 
dt' ds' 



+ / G(r,t\r',U)v^(r',U)dr'. (16) 



The last term can be made zero by choosing an appropri- 
ate initial condition. The first term gives the influence 
of sources on the evolution of the v(x,t) field. The sec- 
ond term incorporates the influence of boundary condi- 
tions. To apply Dirichlet boundary conditions one finds 
a Green's function which is zero at the left boundary. 
Since we have a semi-infinite domain, we use the method 
of images to write down the Green's function 



G(r,t\r',t') = 

e -q 2 (t~t') 

~ v/47r(t - T) 



VMt - 7) 



exp 



(r - r') 2 



exp 



'^^2 



jr + r' + 2x(t')) 
4(t-tf) 



(17) 



where the second term is the image of the first and G = 
at r = —x, the boundary. Proceeding with the calcula- 
tion of source effects in Eq. I|16|l . we have, 



v (n) (r,t) = 



(t')Jo VMM 



exp 



(r - r'f 



exp 



(r + r' + 2x{t')Y 
4(*-f) 



pW(r',t')dt'dr'. (18) 



The Green's function terms in Eq. (|18fl above contain 
exponentials of functions of time t' that possess a max- 
ima at some time t' a . Hence most of the contribution to 
the integral comes around this maximum value, which is 
approximately given by t — t' a = \r — r'\/2q for the first 
integral, and t — t' = \r + r' + 2x(t' )\/2q for the sec- 
ond. If the width of the maxima peak is less than these 
time differences, we can take the limit t — > 00 in the inte- 
grals above. Further, if all the source terms are smoothly 
varying functions of t', one could perform a steepest de- 
scent approximation to the integrals above assuming that 
the maxima peak is nearly a gaussian and sharp enough. 
Physically, as a consequence of both reaction and dif- 
fusion, the configuration of fields at present time t in 
Eq. I|14f> is determined by the time behavior of sources 
in an earlier small time window, in which reaction and 
diffusion mechanisms combine to produce the maximum 
rate of change of the v(x,t) field. At all other times, 
either reaction or diffusion is individually dominant and 
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not able to produce a combined high rate of change of 
v(x,t). This is unlike pure diffusion, where all the time 
history of sources is required to give the field configu- 
ration at present time. Performing the steepest descent 
calculation, the time portion of the integral in Eq. I|18() 
can be eliminated and it reduces to, 



,(n) 



(r) 



pe 



-(r+r'+2x)f 



V2^a(r + r' + 2x) + b 



where, 



p {n \r')dr', 
(19) 



<J 
P 
f 
a 
b 



-q\r — r 



/2<7, 



(±±V± 2 +q^)/2q 2 , 
l/Ap, 



q 2 p 



xp , 
2x V 



2xp + 1/2. 



The first factor in the bracelets g dominates the evolution 
far away from the boundary, while the second factor is 
responsible for sensing the boundary. One notices that 
x and x, the velocity and acceleration of the front are 
involved in the reduced Green's functions in Eq. (|19fl . If 
x and x are neglected in the expressions above (justifiably 
so since we are close to the front bifurcation), solving 
Eq. H14|) reduces to solving 



0, 



1,2,3.. 



(20) 



This is what the authors do in Ref. [Tj], although no 
boundary influence is considered. Solving Eq. (|19fl fur- 
ther for n — 1,2,3, and requiring the smoothness of the 
front at all orders, we get the equation for the evolution 



of 



v f = - — 



da dx 
dx dt 

iq 2 aV2 



da dx da dx 
dx dt dx dt 
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where, 



3/V2: 
[(2- 



<)/Aq 3 + Si + S 2 



(21) 



and 



Si 



So = 



pe 



-(r+r'+2x)f 



-x V2y/ a{r + r' + 2x) + b 



,qr 



pe 



-(r+r'+2x)f 



2y/a(r + r' + 2x) + b 



dr\ 



dr'. 



<l 



Hence Si, S2 and so a depend on x, x, x. Neglecting x 
and x in Si and 52, a reduces to, 



Aq 3 



2 + (1 + 2qx)e- 2 i x ' 



(22) 



Far away from the boundary influence, one recovers 
a = 2q 3 as expected. Therefore, one of the effects of 
proximity to the boundary is spatial dependence of the 
critical bifurcation parameter r\ c . Due to this form of er, 
where its dependence on velocity and acceleration is ig- 
nored and its spatial derivative falls off sharply with x, 
all its derivatives in Eq. (|21|l can be neglected. Thus we 
have, 

Aq 2 aV2. . a 2 , 

v f = (r) c -ri)Vf - —v f . (23) 

We now derive the effect of boundary conditions on the 
evolution of the front. Taking the derivative of G with 
respect to r' in Eq. (|17fl and substituting it in Eq. I|16|) . 
the boundary contribution is found to be 



v c [x(t')]e 



-« 2 (t-t') 



2n 1 ' 2 {t-t'f /2 



x [r + x(t')] exp 



iff 



A(t-t') 



dt 1 . (24) 



Now we are only interested in the contribution of the 
boundary terms at the front position, 0i| r= o- This ex- 
tra term gets added on to Vf, the value of the v field 
at the front position, thus incorporating the influence of 
the specific Dirichlet boundary condition on the front 
velocity. From now on <pi will stand for 4>i\ r= Q. 

Nonsteady front motion represented by v'f involves 
time derivatives, hence the time derivative of <pi gives 
the influence of the boundary condition in accelerating 
fronts, 



dt 



v c [x(t')]e 



-q 2 (t-t') 



2n 1 ' 2 {t-t'f /2 
3x(t') 



q 2 x{t') 



(t - f) 4(i - t>) 




Incorporating these boundary effects in Eq. vj 
Vf + <pi , and v'f — > v'f + 4>2 , we get 



Vf = <t>2 + 



Aq 2 aV2. 

2 I" ^ (??c 



v)(v f + <f>i) - — (vf + faY 



(26) 



Eq. © and Eq. I|26l) constitute the coupling of the two 
degrees of freedom c and x or equivalently Vf and x in 
the presence of a spatial inhomogeneity introduced by 
Dirichlet boundary conditions. 

One of the stationary points of the system above is 
{x,Vf) = (— ln(l — Vbq 2 )/q,0). This is the point closely 
connected with the dynamics we now describe. We look 
at the parameter Vb in three different regimes, Vb < 0, 
< Vb < 1 /q 2 , and Vb > l/q 2 - The region of physical 
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0.015 




012345678 
X 

FIG. 5: The dark curves are the nullclines, grey curves are the 
solutions to Eq. (1261 with different initial conditions. These 
are the typical flows in the regime Vt, < 0, where incoming 
fronts always bounce. There are no fixed points and vt = 
-0.1, ol = 9.9, e = .001, 8 = 1.0. 
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FIG. 6: Flows in the coexistence region < vt < 1/q 2 . Dark 
curves are the nullclines, which intersect at the fixed point. 
Grey curves are the solutions to Eq. (1261 , with different initial 
conditions. The dashed lines represent the invariant manifolds 
separating basins of attraction and bouncing, vt = 0.088, 
ol = 9.9, e = .001, S = 1.0. 



interest is a; > 0. Vf > implies a front moving towards 
the boundary and Vf < is a front moving away. 

Although, the convolution integrals <fi± and 02 repre- 
sent the complete influence of boundary conditions, we 
wish to approximate them to obtain a simpler picture 
that preserves the qualitative features of the complete 
integrals. The integrals involve exponentials of functions 
that have a maxima at time t' a , and hence most of the 
contribution is around this maxima. For small front ve- 
locities this maxima is given by t — t' a — x/2q. The width 
of this maxima peak is given by (x/2q)2. If x/2q >> 
(x/2g)5, which it is for a very sharp peak, a steepest 
descent approximation can be made. A greater inequal- 
ity implies a better approximation. The approximation 
gives, 4>\ — —v c (x)e~ qx and (p2 — d<j>i/dt = —3(p\q/x, 
where it is again assumed that the velocities are small. 

In the regime Vb < the fixed point is in the negative 
x region and is an unstable spiral. The term involving 
1 jx in the expression for (j> 2 prevents flows from crossing 
the negative x region to the positive one and vice versa. 
Therefore the fixed point does not influence the x > 
flows. Figure shows the nullclines and the typical flows. 
The grey flow curves show the turning around of fronts 
at the boundary. It is notable that the nullclines are 
not followed well by the flow curves and cannot predict 
the dynamics of front reversal. Generally, the flows will 
agree better with the nullclines when the relaxation rate, 
determined by rj c — r], is larger, although the jump from 
one nullcline to another can only be explained by the 
dynamical equations. 

As one increases Vb and enters the < Vb < 1/q 2 
regime, the fixed point crosses over into the x > re- 
gion. The fixed point now influences the flows close to 
the boundary. Instead of being an unstable spiral it now 
is a saddle with two distinct real eigenvalues, giving rise 



to stable and unstable manifolds. The eigenvalues are 



A D ± y/Xl + 4X X 



where X Q = Aq 2 a^/2/9{-q c — 77), and 



(1 - v b q 2 ) 



V2~V 



In (1 - v b q 2 ) 



(27) 



(28) 



As Vb —> 1/q 2 , the fixed point moves away from the 
boundary towards positive infinity and A+ — > A , which 
is the eigenvalue for an unstable Ising wall far from the 
boundary influence. Also, in the same limit, A_ — > 0, 
where the zero eigenvalue is associated with spatial ho- 
mogeneity (translational invariance). This explains the 
critical slowing down observed in the last section. Trajec- 
tories wandering close to this fixed point near criticality 
(vb — * 1/q 2 ) will rebound or trap on a slower time scale, 
compared to a relatively faster dynamics when the fixed 
point is further away from criticality. 

The time scales close to the fixed point are controlled 
by A, which has two constituents, A and X x . X is asso- 
ciated with the slow time scale T = c 2 t, which depends 
on the distance to the front bifurcation ry c — 77, and can 
be made arbitrarily small. Therefore, close to the front 
bifurcation, A Q can be neglected in Eq. (|27() and the eigen- 
values reduce to |A±| = %/AT. This is the new time scale 
determined solely by the influence of boundary condi- 
tions and is the dominant time scale in the nonadiabatic 
limit of extremely slow velocities. Bloch wall trajectories 
close to the saddle, which either trap or bounce, evolve 
on this time scale. 

Typical flows are plotted in Fig. EI The triangle shows 
the fixed point. The dashed lines are the invariant sets, 
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FIG. 7: Flows in the regime vt, > 1/q 2 . Again there are 
no fixed points present and all the flows get trapped at the 
boundary. v b = 0.1, al = 9.9, e = .001, 5 = 1.0 

0.015 i j i i 1 1 1 1 1 1 1 
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X 

FIG. 8: This shows the actual simulation of Eq. Q for Vb = 
0.088, al = 9.9, e = .001, 5 = 1.0 The thick dark lines are the 
nullclines. The thin lines are the trajectories. 

with arrows showing the direction of flow. The invari- 
ant sets separate basins of attraction of flows towards 
the boundary and basins of reflection away from it. This 
explains the coexistence region. If the initial velocity 
and position of the front is inside the attraction basin, it 
gets trapped at the boundary, if not, it rebounds. The 
numerically obtained transition curve in the last section 
is a manifestation of this behavior, where initial condi- 
tions are fixed and varying the boundary values leads to 



a crossover of the initial condition from a basin of at- 
traction to that of repulsion. It should be noted that 
the coexistence region can only be explained by the pres- 
ence of the fixed point and the dynamics associated with 
it. An adiabatic analysis relying on nullclines is not the 
complete picture. 

For Vb > 1/q 2 the fixed point no longer exists. Figure[7| 
shows the flows in this regime. Incoming fronts always 
get trapped. Since no fixed points are present, the flow 
qualitatively does what the nullclines do, as is the case 
in the v & < regime. 

Summarizing, transition from bouncing to trapped 
fronts is governed by a fixed point close to the bound- 
ary. This fixed point gives rise to the coexistence be- 
havior, and is absent in regimes where only trapping or 
bouncing occurs. We conclude our analysis by pointing 
out that solutions of Eq. I)23|l , with approximated 4>i and 
4>2, agree well qualitatively with the solution of the FN 
model Eq. |@} plotted in Fig. |S] 



V. CONCLUSION 

We have studied Bloch front motion in the CGL and 
FN models in the presence of spatial inhomogeneity intro- 
duced by Dirichlet boundary conditions at the boundary. 
We have shown similar features in the behavior of Bloch 
fronts close to the boundary in both systems. There is 
a transition from trapping (annihilation) to bouncing of 
incoming Bloch fronts for both systems as a function of 
Dirichlet boundary values. Also for certain boundary 
values trapped and bouncing Bloch fronts coexist. 

In the sharp front and slow velocity regime for the FN 
model, we were able to give a mathematical mechanism 
for the trapping, bounce and coexistence of the two. This 
involves the coupling of the front velocity c and front 
position x close to the boundary. Essentially Dirichlet 
boundary conditions act as a barrier to Bloch fronts com- 
ing in. It is shown that this barrier may or may not be 
penetrated depending upon how fast or far away the in- 
coming Bloch wall is created. It is never penetrated, if 
Vb < 0. It is always penetrated if Vb > 1/q 2 - 

One could also derive how c and x are coupled close to 
a boundary for zero flux boundary conditions, by using a 
modified Green's function, with a zero derivative at the 
boundary. Zero flux boundary conditions 4] show rich 
behavior too, (breathing fronts) and the dynamics can 
be explained by such a calculation. 
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